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Abstract 

The telegraph process X(t), t > 0, (Goldstein, 1951) and the geometric telegraph 
process S(t) = sq exp{(/Lf— ^a 2 )t+aX(t)} with [i a known constant and a > a parame- 
ter are supposed to be observed at n+1 equidistant time points U = iA n , i = 0, 1, . . . , n. 
For both models A, the underlying rate of the Poisson process, is a parameter to be 
estimated. In the geometric case, also a > has to be estimated. We propose differ- 
ent estimators of the parameters and we investigate their performance under the high 
frequency asymptotics, i.e. A n — ► 0, nA = T < oo as n — > oo, with T > fixed. The 
process X(t) in non markovian, non stationary and not ergodic thus we use approxima- 
tion arguments to derive estimators. Given the complexity of the equations involved 
only estimators on the first model can be studied analytically. Therefore, we run an 
extensive Monte Carlo analysis to study the performance of the proposed estimators 
also for small sample size n. 

key words: telegraph process, discretely observed process, inference for stochastic 
processes. 

MSC: primary 60K99; secondary 62M99 



1 Introduction 

The random motions with finite velocity represent an alternative to diffusion models defined 
by means of stochastic differential equations. The prototype of these models is the telegraph 
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process (see Goldstein, 1951 and Kac, 1974) that describes the position of a particle moving 
on the real line, alternatively with constant velocity +v or —v. The changes of direction 
are governed by an homogeneous Poisson process with rate A > 0. The telegraph process 1 is 
defined by 

X(t) = V(0) [ (-l) N ^ds, t>0, (1.1) 



where the initial velocity V(0) assumes the values ±u with equal probability and {N(t),t > 
0} is the Poisson process on [0, t]. We consider a particle initially located on the real line at 
point xq. 

The explicit conditional density function of the process X(t) has been obtained by Ors- 
ingher (1990) and reads 

p{x,t; x Q ,0) (1.2) 
e- xt f f X r^-„ ; nr\ d , fX 



2v \ \v J dt 

e -\t 

H — {8{x - X - Vt) + 8{x - Xq + Vt)} 



XI -\JV 2 t 2 - (x - Xq) 2 + —Iq \-^/v 2 t 2 - (X - Xq) 2 > X{|<D-*o|<«t} 



-Xt 



- A \ vXth (^v 2 t 2 - (x-x ) 2 ) 

A/ ( ~Vv 2 t 2 - (X-X ) 2 + 9j9 === \ X{\x-x \<vt} 



2« " / y/v 2 t 2 -(x- Xq) 



e -xt 

H — {5(X — Xq — Vt) + 5(X - Xq + Vt)} 



2v 



^ [ f\ \ vXth {^u t (x,x ) I I e -xt 

XI I -^/u t (x,x ) I + } X{u t (x,x )>o} + -j-$(ut( x > x o)) 

\ V J \/U t {X,Xo) 



where {x : \x — xq\ < vt}, u t (x, Xq) = v 2 1 2 — (x — Xq) 2 and 



w = 



x \2k+v 1 



; (i 2 J k\T(k + v + l) 



\x\ < co, are; a; < 7r, 



is the modified Bessel function with imaginary argument. Note that the second term in 
equation ()1.2|) represents the singular component of the distribution of (jl.lj) . of the position 
of the particle at time t. Indeed, if no Poisson events occur in the interval [0,t], we have 
that P {X(t) = +vt} = P {X(t) = -vt} = \e~ xt . 

Many authors analyzed over the years the telegraph process, see for example Orsingher 
(1985, 1990), Foong and Kanno (1994), Stadje and Zacks (2004). Di Crescenzo and Pellerey 
(2002) proposed the geometric telegraph process as a model to describe the dynamics of the 
price of risky assets S(t). In the Black-Scholes (1973) - Merton (1973) model the process 
S(t) is described by means of geometric Brownian motion 

S(t) = s exp{at + aW(t)}, t > 0. (1.3) 



1 In the literature, this process is alternatively called the telegraph process or the telegrapher's process. 
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where W(t) is a standard Brownian motion and a = fi — |cr 2 . Di Crescenzo and Pellerey 
(2002) assume that S(t) evolves in time according to the following process 

S(t) = s exp{at + aX(t)}, t > 0. (1.4) 

where X(t) is the telegraph process. Given that X(t) has bounded variation, so is S(t) in 
equation p. 4jl . This seems a realistic way to model paths of assets in the financial markets. 
Mazza and Rulliere (2004) linked the process (jl.l)) and the ruin processes in the context of 
risk theory. Di Masi et al (1994) propose to model the volatility of financial markets in terms 
of the telegraph process. Ratanov (2004, 2005) propose to model financial markets using a 
telegraph process with two intensities \± and two velocities c±. While such markets may 
admit an arbitrage opportunity, linking opportunity velocities and interest rates, the author 
proves that the market becomes arbitrage-free and complete. An analogous of the Black and 
Scholes equation is established as well. 

The aim of this paper is the estimation of the parameter A when {X(t),0 <t< T} is 
observed at equidistant times = t < ... < t n (and also a for discrete observations from 
the process (|1.4)l ). We assume that U = iA n , i — 0, . . . , n, hence nA n = T. The asymptotic 
framework is the following: nA n = T fixed and A n — ► as n —>■ oo. Sometimes we will 
use A instead of A n to simplify the writing. When the telegraph process X(t) is observed 
continuously then N(T)/T is the optimal estimator of the parameter A as this statistical 
problem will be equivalent to the one of the observation of the whole Poisson process 2 
on [0, T\. This situation corresponds indeed to the limiting experiment in our asymptotic 
framework. This asymptotic framework is usually known as the "high frequency" scheme in 
the literature on estimation from discrete time observations of processes solution to stochastic 
differential equations of the form 

dY{t) = b(Y(t), 6)dt + a(Y(t), 9)dW(t), (1.5) 

This field has been an active research area during the last twenty years. The reader can con- 
sult Sorensen (2004) for a review on estimation techniques recently appeared in the literature 
since the seminal papers of Le Breton (1976) and Florens-Zmirou (1989): e.g. estimating 
functions, analytical and numerical approximations of the likelihood function, MCMC meth- 
ods, indirect inference, etc. Unfortunately such methods are not directly applicable in our 
case because the telegraph process is not ergodic or stationary nor Markovian. The main 
idea in the paper is to consider the observed increments of the process X(iA) — X((i — 1) A) 
as n copies of the telegraph process up to time A and treat them as if they were independent 
(which is untrue). From this idea we build an approximated likelihood and score function 
from which we derive estimators. We further propose least squares estimators. Equations 
emerging in connection with the telegraph process are always complicated to treat and closed 
form results are quite rare in the literature. This also happens in our work, thus in some 
cases we rely on numerical simulations to study the properties of the estimators. 

It is worth mention that, up to our knowledge, the only references about estimation 
problems for the telegrapher's processes are Yao (1985) and Iacus (2001). The first author 
considers a the problem of state estimation of the telegrapher's process under white noise 

2 For more details on parametric estimation for Poisson process see Kutoyants (1998). 
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perturbation and studies performance of nonlinear filters. The second paper is about the 
estimation of the parameter 9 of the non-constant rate Ag(t) from continuous observations 
of the process. 

The paper is organized as follows. In Section El we will introduce the approximating 
likelihood function 

LnW= n^lPP^n^-lA-l) (1-6) 

where p(X iy A n ; X^, t^i) is defined by (jl.2|) . We will then study the asymptotic properties 
of estimator 

X n = argmaxL n (A) 

A>0 

(which is interpreted as an approximated maximum likelihood estimator) and the properties 
of the estimator obtained on the approximated score function whose properties are easier to 
study, i.e. the quantity A n satisfying 

= 

A=A n 

We also propose a least squares estimator for A based on the second moment of the pro- 
cess X(t). Inference problems and estimators for the parameters A and o of the geometric 
telegrapher's process are considered in Section El Finally, section 0] contains a Monte Carlo 
analysis to study empirically the behavior of the estimators in a finite sample context (i.e. 
non asymptotically). 



d 

— log L n (X) 



2 The scheme of observation and the asymptotics 

We assume that the telegraph process {X(t),0 <t< T}, with X(0) = x = 0, is observed 
only at discrete times < t\ < • • • < t n = T, with ti = iA n , i = 0, . . . , n hence nA n = T. We 
use the following notation to simplify the formulas: X(ti) = X(iA n ) = Xj. The asymptotic 
is considered as n tends to infinity under the conditions A n — > and nA n = T. The interest 
is in the estimation of the parameter A whilst v is assumed to be known. 

As mentioned in the Introduction, if one can observe the whole trajectory, A can be 
estimated as N(T)/T where N(T) is the number of times the process switches its velocity 
during the interval [0, T] which is, of course, the number of Poisson events counted in [0, T]. 
This is certainly the best estimator of A and it is indeed our target. 

The estimation of v is always an uninteresting problem as, if there are no switchings 
in [(i — 1)A„, iA n ] then X^ — X^i = vA n , hence if A n is sufficiently small, there is high 
probability of observing N(t i+ i) — N(ti) = then v can be estimated (actually calculated) 
without error. 

The process X(t) itself is not Markovian. On the contrary, the two dimensional process 
(X(t),V(t)) has the Markov property but a scheme of observation in which one is able to 
observe both the position and the velocity of the process at discrete time instants is not 
admissible, so we can rely only on the observation of the X(t) component. Hence we cannot 
write an explicit likelihood of the process in the form of a product of transition densities as, 
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for example, for the case of diffusion processes. Another unfortunate fact about the telegraph 
process is that it is not even stationary (at the second order) as 

EX(t) = (2.1) 

and 

^w-tO" 1 ^) < 2 - 2 > 

(see Orsingher (1990)) nor it posseses an ergodic property, so we cannot use the same ap- 
proach as proposed, e.g., in Sorensen (2000). We then need an approximation argument as 
follows. 



2.1 An approximation of the likelihood function 

By taking into account the distribution (jl.2j) . we approximate the likelihood of the process 
with the following function 

n 

L n (X) = L n (\\X ,X 1 ,...,X n ) = Y[p(X i ,A n ',X i -. 1 ,U-i) (2.3) 

i=i 

= n{— W^) + — p(»->°) 



1 = 1 



e -AA„ 



8(u n>i = 0) 



where u n>i = u n {X, h Xi_i) = v 2 A 2 n - (Xi - X^) 2 . 

The density p(Xi, A n ; Xi-i, U-i) appearing in (|2.H|) is the probability law of a telegraph 
process initially located in Xj_i, that reaches the position Xi at time £j. To build our 
approximation of the likelihood, we consider the observed increments X$ — X^x as n copies 
of the process X(A) (i.e. the process X(t) up to time A) and we treat them as if they were 
independent. This is of course untrue, but results below show that our is not a bad idea. It 
is clear that (J2.3)) is equivalent to 

Ln{x) = {—) n— |M^J + — ^=r^j (2 - 4) 

e -XnA n l jt- J fX \ vXAJl (±VU^) 
= ~^Tn ^11 ^0 + 1= 

where n + is equal to the number of approximating telegraph processes with at least one 
change of direction (see Figure 1). 

/ _ AAn \ n-n+ 

In the expression (|2.4|1 . the factor ( e J concerns the singular part of the densities 

p(Xi, A n ; while the product represents the absolutely continuous components 

of the distributions of the telegraph processes. Note that for increasing values of A, the 
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Figure 1: Discrete time sampling of the telegraph process. Between X3 and X4 no Poisson 
event occurred with probability e AA " , hence the singular part of the approximated likelihood 
emerges in the L n (X). For this example trajectory n = 5 and n + = 4. 

absolutely continuous component of (|2.4j) has a bigger weight than the discrete component; 
viceversa for small values of A this has consequences in the performance of the estimators as 
shown in Section |3J Figure Q shows how the two components of the function L n (X) emerge 
for this scheme of observation. 

2.2 Estimators on the approximated likelihood 

One estimator that can be derived from L n (X) is the following quantity 

X n = argmaxL n (A) (2.5) 

A>0 

Of course, A n is not a true maximum likelihood estimator as L n (X) is not itself a true likeli- 
hood, nevertheless such estimators can be effective like in Kessler (2000) but unfortunately, 
our model is not ergodic or stationary and the likelihood looks quite difficult to handle. Nev- 
ertheless we are able to prove uniqueness of the estimator. On the contrary, it is not clear 
how to prove asymptotic properties of the estimator directly on ()2.5|) as the term log(2 n t> n+ ) 
diverges in 

logL n (A) = -XnA n - \og(2 n v n+ ) + J^log <^ A/ + ^ \ (2-6) 

i=l I \ V ' \J u n,i I 
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On the contrary the following estimating function 

F(\;X 1 ,...,X n ) = -^logL n (\) (2.7) 

is easier to study. The estimator in this case is the value A n solution of F(X) = 0, i.e. 

A n :F(A = A n ;X 1 ,...,X n ) = (2.8) 

and once again, we borrow the approach from the literature on estimating functions for 
discretely observed diffusion processes. The function ()2.7|) is indeed more tractable and the 
asymptotic properties of the estimator A n can be studied as in the next theorem. We just 
note that, given the uniqueness of A n (proved in Theorem 12 . 21 below) . A n coincides with A n 
hence they share the same properties and they will be treated as one in Table d 

Theorem 2.1 Under the condition nA n = T as A n — > we have that 

An ^Xoo = ■ ( 2 - 9 ) 

Proof. In order to prove (J2.9)) we recall some properties of Bessel's functions (see e.g. §5.7, 
Lebedev, 1972) 

^J n (x) = ^ ^I n -i(x) + J„+i(x) 
and 

iAk ■ u) k 

lim = — , lim Io(k ■ u) = 1, lim iAk ■ u) = 0, lim ^(k ■ u) = 0, 

u^O u 2 u^O u^O u— >0 

Direct differentiation of ()2.(ij) gives 

_ A h (jkgg (i + + (^g + jfe) h gygg + 



In the limit as n — > oo we have, A n — > 0, nA„ — > T, u = ii n> i = t> 2 A 2 — (Xj — Xj_i) 2 — > 
and n + — > N(T). Therefore ()2.10|) converges to 

A 

□ 

In light of Theorem 12.11 we can say that in the high frequency observation scheme, the 
estimator A n tends to the best estimator of the intensity A of an homogeneous Poisson 
process. 

Theorem 2.2 The estimator A n in (|2.5jl is unique. 
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Proof. To prove the result, we show that the second partial derivative of logL n (A) is 
negative. In view of the following property of Bessels functions 



we have that 



I v „ 1 (z)-I v+1 (z) = ^I v (z) (2.11) 



AA„ f „ f\ \ f \ \ 1 A n v f \ 



2 1 u I V / Cf » * l « l v "'lift i f . 

therefore equation ()2.10j) can be rewritten as 

d ^ ^n—t; (1 + A n A) 7 (jVgg + u n<i A A 
-logL.(A) = -nA„ + g ^ ( ^ /o ( A v/ -).,. A! , /; (A^-)) ( 2 - 12 ) 

Differentiating 1)2.12)1 again with respect to A and posing x = ^y/u n< i to simplify the equa- 
tions 3 , we obtain 

d 2 n+ 

—— log L n (A) = V] -2 

9X2 ~[2v 2 \ 2 (^JIo(x) + v A n h(x)) 2 

where the generic term gi is the sum of the following terms 

-2^I-v 3 A n I 1 {x)I {x) < (2.13) 

u n ,i {u nil A 2 - v 2 (2 + A n A + A 2 A 2 ) } I 2 (x) (2.14) 

2u n>i X {-\u nii + v 2 A n (1 + A„A)} If{x) (2.15) 

u n>i A (u n .i X -v 2 A n (1 + A n A)) I 2 (a;)io(aO (2.16) 
Recalling that u n> i = v 2 A 2 l — (Xj — Xj_i) 2 > 0, equation 1)2.14)1 can be rewritten as 

-u nii {v 2 A n \ + (Xi - AVO 2 A 2 } I 2 (x) - 2m„^ 2 I 2 (x) (2.17) 

summation of ()2.15jl and ()2.16)1 gives 

-itn^ {t; 2 A n A + (X,-A l _ 1 ) 2 A 2 } (2 h(xf - I (x)I 2 (x)) (2.18) 

Putting together 



{v 2 A n A + (A, - X^) 2 A 2 } J 2 



and equation ()2.18j) it remains to study the sign of 

I 2 (x) + 2h{xf - I (x)I 2 (x) = I (x)(I (x) - h{x)) + 2h{xf 
which is positive due to the fact that Iq{x) > Iz(x) for positive x (from property 1)2.11)) ). □ 



3 We remind that this x — ^Ju^lX/v is strictly positive because it refers to the n + terms for which u n> i > 0. 



2.3 A least squares estimator 



We already mentioned that the telegrapher's process is not stationary at the second order, 
but we can still think to use (|2.2jl to obtain a least squares estimator on the mean of the 
squared increments in the following way If, as before, we consider the observations X{ — 
as n copies of the process X(A), then looking at ()2.2|) for t — A we have 

v 2 ( 1 - e~ 2XA 
EX 2 (A) = ^ ( A 



A V 2A 

Consider now the sample second moment (or the mean of the square) of the observed incre- 
ments Xi — 



m 2 

n 

i=l 



^(Xt-x^y 



then the following estimator can be considered 

X n = argmm |m 2 - - (A n ^— j | (2.19) 

This estimator will be compared numerically to the estimator A n in Section 



3 Parametric estimation for the geometric telegraph 
process 

Consider the process Y of the observed log-returns 

Yi = log = aA n + a(Xi - X^ x ) 

where Si = S(ti) are discrete observations from the geometric telegraph process (jl.4|) . We 
assume /i to be known, which is usually the case in finance where \i is related to the expected 
return of non risky assets like bonds, etc. The parameters a and A are to be estimated. As 
in the previous sections, we can assume v to be known as well, if not we will show in the 
next paragraph a simple way to obtain it. We assume Yi to be n copies of the process 

Y(A) = aA + aX(A) 

with X(A) =Xi- Xi-i and X(0) = x = 0. Therefore, by (J2HJ), we have 

EY(A) = aA 

and by ()2.2|) we obtain 

v 2 { 1 - p~ 2AA \ 
VarF (A) = (7 2 VarX(A) = a 2 — ( A - J (3.1) 
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A good estimator of the volatility a can be derived from the sample mean of the log returns. 
Indeed, 



1 

Y n = -YY i = aA + -X n (3.2) 
i=i 

and 



EY n = «A+ ^EX„ =aA = - ^<j 2 j A (3.3) 

again by (|2.1|) and for the properties of the log-returns. From ()3.3|) we have that 

2 f VY n 
a = 2^- — 

from which the following unbiased estimator of a 2 can be derived 

*»= 2 ("-!)■ 

Therefore, a reasonable estimator of cr is 



*» = a/2(a*-^) (3.4) 

which not always exists because there is no guarantee that fx > Y n /A. Moreover, it should 
be noticed that in practice, given /i, a and A the estimator essentially depends on the last 
value of the telegraph process X n . In fact, we can write ()3.4|) in terms of the telegraph 
process 



therefore the estimate of a does not depend on n. This is why Table H] reports the same 
value of the estimates for different sample sizes. We then use a n to estimate A making use 
of flUU). Let 



-Y,{Yi-Y n 

then the proposed estimator of A is 



n 

i=l 



X n = arg mm (4 - ft- (A - j j (3.5) 

3.1 Filtering of the geometric telegraph process 

If the velocity v is not known one can proceed as follows: set 

Y - EY 

Z t = 1 n = X l -X l _ 1 = X{A) 
a 
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hence an estimator of the increments of the telegraph process is 



Z l = Y t -^ = X(A), z = l,...,n 

then 

Z\ = Xi, Z2 + Zi = X2, Z3 + Z2 + Z\ = A3, . . . 

where Xi are the estimated states of the underlying telegrapher's process. From these esti- 
mates, one can proceed as in previous sections and estimate both A and v. 

4 Monte Carlo analysis 

To assess the properties of the estimator ()2.5|) and (|2.8J1 for fixed n < 00 we run extensive 
Monte Carlo analysis. Given that numerically A n and A n coincides, we treat them as one 
in the tables. We simulate 10000 trajectories of the telegrapher's process on the interval 
[0,T], T = 500, for different values of A and with v = 1 fixed. Each trajectory has then 
been resampled on a regular grid of n = 50, 100, 500 and 1000 points and the corresponding 
observations have been used to estimate the unknown parameter. The results have been 
collected in Table [T] It emerges that, as expected, the bias tends asymptotically to zero as 
well as the mean square error. Furthermore, bias and variance are strictly correlated to the 
value of the unknown parameter A. This is expected as well because, for fixed n, as the more 
A increases the more Poisson events remain hidden to the observer. For the same experiment 
(and on the same sample trajectories), Tabled reports the performance on the least squares 
estimator A n from equation (|2.19j) . 

Tables E] and |U reports estimates results on the telegrapher's process respectively for the 
estimation of A and a. The paths of the geometric telegrapher's process have been generated 
from the ones of the telegrapher process of Tables Q and El As it can be seen, the estimator 
A n in Table El strongly depends on the quality of the estimate a n (reported in Table HJ). For 
low values of A, in some cases the condition for the existence of a n , i.e. \i > Y n /A, has 
not been fulfilled hence we report the percentage of valid paths over the 10000 simulated. 
For A > 0.75, it seems that A ra performs quite similarly to A ra in terms of bias. This seems 
consistent with the definition of the estimators and the performance of the estimator a n . 
Tables also report the column yMSE. Values under this column are calculated, e.g., as 
follows 



MSE(A n ) 



\ 



1 ™ 



where N is the number of Monte Carlo simulations (N = 10000 in our case) and n is the 
fixed sample size. 

5 Conclusions 

Despite the analytical problems in dealing with the telegraph process, in this paper we 
have shown that statistical inference on the process can be attempted. Moreover, numerical 
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results seem to suggest that also for the geometric telegraph process this attempt might be 
successful which encourages the study of the analytical properties of this process in view of 
financial applications. In particular, a detailed description of the law of the increments of 
the process of the log-returns might be of interest because, numerical evidence, show that 
these have heavy tails but their law do not need a high number of parameters like other 
distributions proposed in the literature (see e.g. Eberlain and Keller, 1995). Moreover, the 
parameters have a direct interpretation as in the standard geometric Brownian motion. This 
will be a topic for future research. 
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A 


Bias 


a/MSE(A) 


min A 


max A 


n 


0.10 


-0 


002 


0.018 


0.04 


0.20 


50 




-0 


001 


0.016 


0.05 


0.16 


100 




-0 


000 


0.014 


0.06 


0.16 


500 




-0 


000 


0.014 


0.05 


0.15 


1000 


0.25 


-0 


on 


0.041 


0.13 


0.47 


50 




-0 


003 


0.031 


0.16 


0.41 


100 




-0 


000 


0.023 


0.16 


0.34 


500 




-0 


000 


0.023 


0.16 


0.35 


1000 


0.50 


-0 


062 


0.092 


0.26 


0.78 


50 




-0 


Oil 


0.059 


0.32 


0.85 


100 




-0 


001 


0.035 


0.37 


0.65 


500 




-0 


000 


0.033 


0.37 


0.63 


1000 


0.75 


-0 


151 


0.175 


0.36 


1.01 


50 




-0 


031 


0.091 


0.47 


1.18 


100 




-0 


001 


0.048 


0.58 


0.96 


500 




-0 


000 


0.043 


0.60 


0.92 


1000 


1.00 


-0 


264 


0.283 


0.45 


1.23 


50 




-0 


064 


0.128 


0.62 


1.53 


100 




— 


r\r\ -i 

001 


0.058 


0.79 


1.26 


500 




-0 


001 




n si 


1 99 


i nnn 

1UUU 


1.50 


-0 


546 


0.558 


0.58 


1.48 


50 




-0 


162 


0.227 


0.90 


2.16 


100 




-0 


003 


0.080 


1.24 


1.86 


500 




-0 


001 


0.066 


1.27 


1.77 


1000 


2.00 


-0 


874 


0.882 


0.75 


1.65 


50 




-0 


298 


0.357 


1.11 


2.67 


100 




-0 


006 


0.106 


1.63 


2.47 


500 




-0 


000 


0.083 


1.68 


2.33 


1000 



Table 1: Empirical performance of the estimator A n of ()2.5|) for different values of the 
parameter A and different sample size. The time horizon T is fixed to 500. Results over 
10000 Monte Carlo paths of the telegraph process. See text for more details. 
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A 


Bias 


a/MSE(A) 


min A 


max A 


n 


0.10 


0.002 


0.022 


0.04 


0.24 


50 




0.001 


0.018 


0.05 


0.18 


100 




0.000 


0.016 


0.05 


0.16 


500 




-0.000 


0.016 


0.05 


0.16 


1000 


0.25 


0.007 


0.051 


0.13 


0.55 


50 




0.003 


0.037 


0.14 


0.41 


100 




0.000 


0.026 


0.16 


0.36 


500 




0.000 


0.025 


0.17 


0.35 


1000 


0.50 


0.018 


0.106 


0.25 


1.13 


50 




0.007 


0.070 


0.32 


0.93 


100 




0.001 


0.040 


0.35 


0.66 


500 




0.000 


0.037 


0.35 


0.65 


1000 


0.75 


0.028 


0.161 


0.40 


1.72 


50 




0.012 


0.105 


0.45 


1.32 


100 




0.001 


0.054 


0.56 


1.01 


500 




0.000 


0.048 


0.58 


0.96 


1000 


1.00 


0.040 


0.219 


0.53 


2.42 


50 




0.017 


0.141 


0.62 


1.80 


100 




0.002 


0.066 


0.73 


1.27 


500 




n nm 

U.UUl 


U.UO ( 


n 77 


l.ZO 


1 nnn 

1UUU 


1.50 


0.059 


0.329 


0.72 


3.00 


50 




0.028 


0.218 


0.92 


2.78 


100 




0.003 


0.093 


1.19 


1.89 


500 




0.001 


0.075 


1.25 


1.79 


1000 


2.00 


0.080 


0.412 


1.05 


3.00 


50 




0.035 


0.290 


1.22 


3.00 


100 




0.006 


0.120 


1.59 


2.46 


500 




0.003 


0.095 


1.66 


2.39 


1000 



Table 2: Empirical performance of the estimator A n of ()2.19j) for different values of the 
parameter A and different sample size. The time horizon T is fixed to 500. Results over 
10000 Monte Carlo paths of the telegraph process. See text for more details. 
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A 


Bias 


X /MSE(A) 

V \ / 


min A 


max A 


% valid 


n 


0.10 


0.018 


0.107 


0.00 


1.01 


96 


50 




0.042 


0.157 


0.00 


1.12 


96 


100 




0.301 


0.609 


0.00 


3.40 


96 


500 




0.634 


1.189 


0.00 


6.48 


96 


1000 


0.25 


0.009 


0.132 


0.00 


0.91 


99 


50 




0.006 


0.161 


0.00 


0.95 


99 


100 




0.123 


0.427 


0.00 


2.27 


99 


500 




0.320 


0.776 


0.00 


4.00 


99 


1000 


0.50 


0.022 


0.185 


0.00 


1.64 


100 


50 




0.010 


0.186 


0.00 


1.27 


100 


100 




0.031 


0.395 


0.00 


2.03 


100 


500 




0.139 


0.635 


0.00 


3.33 


100 


1000 


0.75 


0.031 


0.239 


0.08 


1.99 


100 


50 




0.014 


0.215 


0.00 


1.90 


100 


100 




0.002 


0.394 


0.00 


2.25 


100 


500 




0.049 


0.612 


0.00 


3.28 


100 


1000 


1.00 


0.045 


0.304 


0.31 


2.68 


100 


50 




0.021 


0.256 


0.26 


2.24 


100 


100 




-0.000 


0.388 


0.00 


2.65 


100 


500 




0.012 


0.608 


0.00 


3.69 


100 


1000 


1.50 


0.061 


0.419 


0.57 


4.22 


100 


50 




0.028 


0.330 


0.59 


3.22 


100 


100 




-0.004 


0.380 


0.03 


3.00 


100 


500 




-0.014 


0.580 


0.00 


3.97 


100 


1000 


2.00 


0.094 


0.529 


0.80 


6.49 


100 


50 




0.038 


0.402 


0.93 


4.16 


100 


100 




0.005 


0.380 


0.64 


3.46 


100 


500 




-0.002 


0.536 


0.00 


4.08 


100 


1000 



Table 3: Empirical performance of the estimator A n of ()3.5j) given the estimate a n (see Table 
HJ, for different values of the parameter A and different sample size. The time horizon T is 
fixed to 500. Results over 10000 Monte Carlo paths of the geometric telegraph process, val 
= % of valid cases, i.e. simulated paths such that \i > Y n /A. See text for more details. 



16 



A 


Bias 


V / MSE(a) 


min a n 


max(X n 


% valid 

cases 


n 


0.10 


-0.008 


0.138 


0.01 


0.85 


96 


50, 100, 500, 1000 


0.25 


-0.007 


0.094 


0.02 


0.74 


99 


50, 100, 500, 1000 


0.50 


-0.003 


0.065 


0.09 


0.69 


100 


50, 100, 500, 1000 


0.75 


-0.002 


0.053 


0.22 


0.66 


100 


50, 100, 500, 1000 


1.00 


-0.002 


0.045 


0.30 


0.66 


100 


50, 100, 500, 1000 


1.50 


-0.002 


0.037 


0.33 


0.64 


100 


50, 100, 500, 1000 


2.00 


-0.001 


0.031 


0.36 


0.61 


100 


50, 100, 500, 1000 



Table 4: Empirical performance of the estimator a n of ()3.4|) for different values of the param- 
eter A, different sample size and o = 0.5. The time horizon T is fixed to 500. Results over 
10000 Monte Carlo paths of the geometric telegraph process. '% valid cases' = percentage 
of valid cases, i.e. simulated paths such that the estimator of o exists, i.e. fi > Y n /A. See 
text for more details. 



17 



